SHARP HEAT KERNEL ESTIMATES IN THE FOURIER-BESSEL SETTING 
FOR A CONTINUOUS RANGE OF THE TYPE PARAMETER 
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Abstract. The heat kernel in the setting of classical Fourier-Bessel expansions is defined by an 
oscillatory series which cannot be computed explicitly. We prove qualitatively sharp estimates 
of this kernel when the associated type parameter is not less than —1/2. Our method relies on 
establishing a connection with a situation of expansions based on Jacobi polynomials and then 
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transferring known sharp bounds for the related Jacobi heat kernel. 



<^ ■ 1. Introduction 



Let J v denote the Bessel function of the first kind and order v > —1, and let {A„ jt , : n > 1} 
be the sequence of successive positive zeros of J v . The Fourier-Bessel heat kernel is given by the 
oscillating sum 
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(1) G t (x, y) = 2{xy) > exp ( - i A ) — — . 

Our main result is the following (see the end of this section for the notation). 

OV 



Theorem A. Assume that v > —1/2. Given any T > 0, we have 
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G?(s,y) ^ [tVxy]— 1/2 



(l-x)(l-y) a 1 
t 



1 / (x-yf 
—p exp — c 



uniformly in x,y G [0,1] and < t < T. Moreover, 

G v t (x,y) -exp(-a^) (l-x)(l-y), 
uniformly in x,y £ [0,1] and t >T. 

Thus we obtain a qualitatively sharp description of the kernel ([T|). The restriction on v 
is imposed by a similar restriction in the Jacobi setting from which the heat kernel bounds are 
transferred. Nevertheless, it is natural to conjecture that the same estimates hold for all v > — 1. 

The main contents of Theorem [X] are of course the short time bounds. Proving them by a 
direct analysis of the heavily oscillating series in ([T]) is practically impossible. Note that the 
Bessel function J v is transcendental in general, and can be expressed by means of elementary 
functions only if the index v is half-integer, i.e. v = k/2 for some integer k. Furthermore, the 
zeros of J v are known explicitly only when v = ±1/2. Note also that the order of magnitude 
of the numbers v appearing in the exponential factor is n 2 as n — > oo; in particular, the 
asymptotic distribution of A 2 v is not linear. 
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The behavior of the Fourier-Bessel heat kernel does not seem to have been studied before, 
except for our previous paper [3]. More precisely, in Theorem 3.3] we derived qualitatively 
sharp estimates for G"(x,y) under the assumption that v is a half- integer not less than —1/2 
(actually, this restriction is essential only for the short time bounds, see Theorem 3.7]). We 
also conjectured analogous estimates for general v > —1, see [31 Conjecture 4.3], and Theorem [XI 
confirms this for all v > —1/2. However, the methods we use in the present paper are much 
different from those applied in [3]. Here the main argument is based on a relation we establish 
between the Fourier-Bessel setting and the framework related to Jacobi 'functions'. This makes 
it possible to transfer qualitatively sharp estimates of the Jacobi heat kernel obtained recently 
by Nowak and Sjogren [3]. 

Our principal motivation for investigating G\ (x, y) comes from an interest in harmonic analysis 
related to Fourier-Bessel expansions, see the references in [31 Section 1]. Additional motivation 
emerges from the probabilistic interpretation of the Fourier-Bessel heat kernel, see [21 Section I]. 
Namely, it is a transition probability density for the time scaled Bessel process X^ 1 ^ 2 on (0, oo) 
(with reflecting barrier at x = when — 1 < v < 0, see [H Appendix I, Section 21]) killed upon 
leaving the interval (0, 1). We also note that for v = d/2 — 1, d > 1, the kernel G"(x, y) provides 
a fundamental solution to the standard heat equation in the c?-dimensional Euclidean unit ball 
with Dirichlet boundary condition and a radial initial condition. 

There are two simple cases occurring already in Fourier's works, at least implicitly. These 

I /2 1/2 

are G t {x,y) and G t {x,y), and the two kernels can be written by means of non-oscillating 
series, see [31 Section 4]. The argument is based on the periodized Gauss- Weierstrass kernel 
and simple boundary-value problems for the classical heat equation in an interval. No more 
elementary representations in these cases seem to be possible, and this suggests that a general 
closed explicit formula for G"(x,y) does not exist. The estimates in Theorem lAl are therefore a 
natural and desirable substitute for an exact expression. 

When writing estimates, we use the notation X < Y to indicate that X < CY with a positive 
constant C independent of significant quantities. We write X ~ Y when simultaneously X < Y 
and Y < X. Further, we will also use the notation X ~~ Yexp(— cZ) to indicate that 

C~ X Y exp(-cyZ) < X < Cyexp(-c 2 Z) 

with positive constants C, c\ and c 2 independent of significant quantities. 

2. Fourier-Bessel and Jacobi settings 
Given v > — 1, consider the Bessel differential operator 

LV _ d 2 2v + 1 d 

dx 2 x dx 

on the interval (0, 1). It is well known that L v is symmetric and nonnegative on C 2 (0, 1) C 
L 2 (dfi u ), where \x v is a power measure in the interval [0, 1] given by 

dpuix) = x 2u+l dx. 

A classical orthonormal basis in L 2 (d[i u ) consisting of eigenfunctions of L u is the Fourier-Bessel 
system {4> u n : n > 1} (cf. [7, Chapter XVIII]) defined by 

= .. x- v J v {\ n , v x), n > 1, x G [0, 1]; 

\Jv+\\An,v)\ 
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here the value <^(0) is understood in the limiting sense. We have 

= A^X, n>l. 

The self-adjoint extension of L v associated with the Fourier-Bessel system is defined naturally 
in L 2 (d/i u ) via its spectral decomposition as 

oo 
n=l 



on the domain 



BomL u = [/ € L\dfi u ) : £ |A^</,K>^J 2 < oo}. 

^ n=i J 



The semigroup generated by — L v has the integral representation 

exp ( - tL v )f(x) = I GUx, y)f(y) d^(y), 



where the kernel is given by (JTJ) . 

By technical reasons, we also need to consider a closely related setting emerging from incor- 
porating the measure [i v into the system {</>^}. In this way we derive the Fourier-Bessel system 
Wn ■ n > 1}, 

M(x) = x v+1 / 2 <%(x), n>l, a; €[0,1], 

which for each v > —1 is an orthonormal basis in L 2 (dx); here and elsewhere dx stands for 
Lebesgue measure in the interval [0, 1]. This system consists of eigenfunctions of the differential 
operator 

d 2 l/A-v 2 

^ ^ dx 2 x 2 ' 

more precisely, we have 



= A^C n > 1. 

The operator L, u is symmetric and nonnegative on C 2 (0, 1) C L 2 (dx). Similarly as in the previous 
setting, there is a natural self-adjoint extension of h u in L 2 (dx), denoted by h u , whose spectral 
decomposition is given by the The associated heat semigroup exp(— £L V ) possesses an integral 
representation (with integration against dx), and the relevant kernel is 

(3) K»(x,y) = (xyy +1 / 2 GUx,y). 

Finally, we introduce the Jacobi setting based on Jacobi trigonometric 'functions', see [U Sec- 
tion 2], scaled to the interval [0, 1]. Let P£ , k = 0, 1, 2, . . ., be the classical Jacobi polynomials 
with type parameters a, ft > —1, as defined in Szego's monograph [6]. Define 

®k ( x ) =c k [ sm ^) l COS irJ k (cosirx), k>0, xe[0, 1J, 

with 

a,p _ fn (2k + a + /3 + l)T{k + a + g + l)r(fc + 1) N 1/2 
T(k + a + l)r(fc + /3 + 1) 
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where for k = and a + (3 = — 1 the product (2k + a + f3 + l)T(/c + a + (3 + 1) must be replaced 



by r(a + /3 + 2). Then the system {$? ,/3 : A; > 0} is an orthonormal basis in L 2 (dx). Moreover, 



each ^ s an eigenfunction of the differential operator 



J' 



a>p _ d 2 vH(l/4-a 2 ) tH(1/4-^) 



dx 2 4sin 2 (7rx/2) 4cos 2 (7rx/2) ' 
and we have 

r ,/^ = a^/W, fc > o, where A^ = Wfc + ^±1V 



k ls -k ^k ' ^ — u > vvncxc "J. 

Thus J"'' 3 has a natural self-adjoint extension in L 2 (dx) given by 

oo 
k=0 

on the domain 

DomJ^ = {/ € L\dx) : |A^</,*^)| 2 < ooj. 

^ k=o > 

The semigroup generated by — S a ^ in l?(dx) has the integral representation 

exp ( - 0^)/(x) = f G^\x, y)f(y) dy, 
Jo 

where the Jacobi heat kernel is defined by the oscillating series 

oo 

G^(x,y) = £exp ( - tkf)$f(x)$f(y). 



k=0 



Recently, qualitatively sharp estimates of (x,y) were obtained in [3], under the restriction 
ct,/3> —1/2. The result below is a consequence of [U Theorem A], the relation between various 
Jacobi heat kernels [H (3)] and a simple scaling argument. 

Theorem 2.1 ([!]). Assume that a, ft > —1/2. Given any T > 0, we have 



(x,y) ~~ 



Q+l/2 



(l-g)(l-y) • 



W2 1 ( (x - y) 2 
: exp — c- 



uniformly in x,y G [0, 1] and < f < T. 

These estimates are crucial for our reasoning proving Theorem lAl 

3. Proof of Theorem [A] 

We first focus on showing the short time bounds, which constitute the main part of Theorem lAl 
Our argument relies on relating the Fourier-Bessel heat kernel K%(x,y) with the Jacobi heat 
kernel G" (x, y) and then transferring the estimates of Theorem 12.11 This will be achieved by 
proving that the generator of the Fourier-Bessel semigroup {exp(— ilL^)} is a 'slight' perturbation 
of the Jacobi generator with the parameters of type a = v and /3 = 1/2. Then the desired relation 
will follow by the so-called Trotter product formula. Note that this transference method works 
for all v > —1, and the only reason for the restriction on v in Theorem |A] is the analogous 
restriction on a in Theorem 12.11 
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A close connection between the generators is suggested by a relation at the level of differential 
operators. Observe that h u and S u ' 1 ^ 2 differ only by a zero order term, 



IS - JT' 1/2 



v 

4 



7T~ 



4 sin 2 



Moreover, the difference function H u (x) is continuous and bounded on [0, 1], with the value at 
x = understood in the limiting sense. Indeed, as easily verified, we have 



\H v {x)\ < \H"{1) 



1 

4""" 



7T 

T ' 



x e [o,i]. 



In the sequel this bound will be used repeatedly, sometimes implicitly. 

An essential point of the whole method is the coincidence of domains stated below. 

Theorem 3.1. For each v > — 1, 

Domf> = Domih 1 / 2 . 
To prove this we will need a preparatory result. 



Lemma 3.2. Xei ^ > —1. Then 



n > 1, fc > 0, 



where IS is the differential operator given by ([2]). 
Proof. We shall use the divergence form of L", 

L "/(x) = -x- I '- 1 / 2 D^+ 1 L>(^- 1 / 2 /(x))) : 



and integrate by parts. Denote 







h u ^{x)^ k ' 1/2 (x)dx. 



Clearly, the integral here converges since and ^j!' 1 ^ 2 are in L 2 (dx). We have 

■ i 



o 



x^+ 1 / 2 D(x- 1 '- 1 / 2 C(x))<I>I ,1/2 (x) '+ / D(x-"- 1 /2^(x))£>(x-'- 1 / 2 ^ 1/2 (x))x 2 -+ 1 dx 



+ 

o Jo 



■x v+1 l 2 D{x- v -V 2 r n {x))®t /2 {xi\ + x 1/+1 / 2 C(x)^(x-^ 1 / 2 $f /2 (x)) 
jT 1 C(^)x-^ 1/2 J D(x 2 ^+ 1 Z)(x-^ 1 / 2 <I>^ 1/2 (x))) dx 



Since X 3 = (C L " $ 



i/,l/2\ 



it remains to ensure that X 



From the definition of &t it is clear that 
<S>fl 2 (x) = 0{x v+1 '% x^0+, 



and 



r fc y^j-^^ J, -r w , < u , > . 1/2 (x) = O ( 1 - x) , X -»• 1 . 

The estimate |0^(x)| < 1 — x, see the proof of O Theorem 3.7], holds uniformly in x G (0,1) 
and implies 

^(x) = C(x^ +1/2 ), x^0+, and <(x) = 0(1 - x), x -> 1". 
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Further, using D(z~ u ' J u (z)) = —z~ v J v+ \{z) (cf. [7, Chapter III, Section 3-2]), we see that 

'0(x), x->0 + 



D{x- u - 1 l 2 r n (x)) = cx~ v J v+l {\ n , v x) 



o(i), s->r 



where c does not depend on x. Moreover, by the differentiation rule for Jacobi polynomials (cf. 
B (4-21.7)]), 



^-P a k \z) = l -(k + a + (3 + l)P^ +1 (z), k > 0, 



(here we use the convention that P^f = 0) it follows that 



D{x- u - 1 / 2 ^ ll2 {x)) = cisin(7rx)P fc l/ + 1 1 ' a/ "(cos(vr2;)) I 



v+1,3/2. 



( sin • 



X 



TTX 



COS ■ 



+ c 2 P, 



",1/2/ 



COS 



(ttx)) 



7T2 7nr 

c 3 P;' 1/2 (cobCttx)) 2 2 



COS 

( : 

cos T? — sin / sin 



s,nfV +1 /2 7TX 

sin T 



.r 



7TX / 7TX 



i/— 1/2 



7TX 



cos ■ 



2 ' 



with ci,C2,C3 independent of x. Therefore 

D(x- v - 1 ^ 1/ \x)) 



0(x), x -> H 
0(1), 

Altogether, these facts imply X\ = X 2 = 0. The lemma follows. 



□ 



Proof of Theorem \3.1[ We demonstrate that DomJ"' 1 / 2 C DomL". Then the other inclusion 
follows by self-adjointness of the two operators. 
Assume that / G DomJT' 1 / 2 and let 

oo 



n=l 



We must show that S is finite. We have 



oo oo 



s = E = E E a** 1 /2 ><*f 

n=l n=l fc=0 

Applying now Lemma 13.21 we can write 

oo oo 2 



5 = E E </> *r ><r ,1/2 + flw /a .# 

n=l £;=0 

OO / CO \ i 

= E ( E ^' 1/2 </> *? /2 )*?'\k) + (E</> *f /2 >^ 1/2 ! 

n=l ^ fc=0 

oo 

= E|<j m/2 /,c> + </^< 



fc=0 



n=l 



This together with Parseval's identity implies 



s<V' 1/2 f\\U dx) +nfHn\h (dx) < 



(dx) 



DC 



and consequently / 6 DomL". 



□ 
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We now recall the Trotter product formula, see Chapter VIII, Section 8]. Let A and B 
be (possibly unbounded) self-adjoint operators on a Hilbert space %. If A + B is essentially 
self-adjoint on Donij4 n DomB, and A and B are bounded from below, then 

(\ m 
exp(-tA/m) exp(-tB/m) ) h, heH, t>0. 
J 

Specifying this general situation to our context we take U = L 2 (dx), A = F' l/2 and choose B 
to be the multiplication operator by H v . Then, in view of Theorem 13. 11 A + B = L". Moreover, 
the Trotter product formula applies and since \H v (x)\ < \H U (1)\, x G [0, 1], it implies 

e^WI exp ( - tF>V 2 )f < exp ( - tl") f < exp ( - tl^ 2 )f, < / G L 2 (dx), 

for each t > 0. Taking now into account that the Jacobi and the Fourier-Bessel heat kernels are 
continuous functions of their arguments (see [H Section 2] , [3l Section 2] ) , we infer that 

e-^^Gf /2 (x,y) < KUx,y) < ^^G^^y), x,y G [0,1], t > 0. 

This combined with Theorem 1 2 . 1 1 and (|3|) justifies the short time bounds stated in Theorem El 
The long time estimates were proved for all v > —1 in [3, Theorem 3.7], under the assumption 
that the threshold T is sufficiently large. However, given any < Tq < T, the already proved 
short time bounds imply 

G%(x,y)~(l-x)(l-y), x,y€[0,l], T <t<T, 

and the desired conclusion follows. 
The proof of Theorem |A] is complete. 

Remark 3.3. A slightly more detailed analysis reveals that for v G [—1/2,1/2] 

e- tHV ^Gf /2 (x,y) < K»(x,y) < e^MG^^y) < G? l/2 (x, y), x,y G [0,1], t > 0, 

and for v £ [-1/2,1/2] 

Gf /2 (x,y) < e- tH ^Gf%,y) < K^y) < e~ tH ^Gf%,y), x,y G [0,1], t > 0. 
Notice that K^ 1 ^ 2 (x,y) = G^ l ^ 2,1 ^ 2 {x,y). 

References 

[1] A. N. Borodin, P. Salminen, Handbook of Brownian motion: facts and formulae, 2nd ed., Birkhauser Verlag, 
Basel-Boston-Berlin, 2002. 

[2] S. Karlin, J. McGregor, Classical diffusion processes and total positivity, J. Math. Anal. Appl. 1 (1960), 
163-183. 

[3] A. Nowak, L. Roncal, On sharp heat and subordinated kernel estimates in the Fourier-Bessel setting, Rocky 
Mountain J. Math., to appear. arXiv: 1111 . 5700vl 

[4] A. Nowak, P. Sjogren, Sharp estimates of the Jacobi heat kernel, preprint 2011. arXiv: 1111 .3145vl 

[5] M. Reed, B. Simon, Methods of modern mathematical physics. I. Functional analysis, Second edition. Aca- 
demic Press, Inc., New York, 1980. 

[6] G. Szego, Orthogonal polynomials, Fourth Edition, Amer. Math. Soc. Colloq. Publ. 23, Amer. Math. Soc, 
Providence, R. I., 1975. 

[7] G. N. Watson, A treatise on the theory of Bessel functions, Cambridge University Press, Cambridge, 1966. 

Instytut Matematyczny, Polska Akademia Nauk, Sniadeckich 8, 00-956 Warszawa, Poland 
E-mail address: adam.nowak@impan.pl 

Departamento de Matematicas y Computacion, Universidad de La Rioja, 26004 Logrono, Spain 
E-mail address: luz.roncal@unirioja.es 



